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Abstract. We consider the reconstruction of optical parameters in a domain of 
interest from photoacoustic data. Photoacoustic tomography (PAT) radiates high 
frequency electromagnetic waves into the domain and measures acoustic signals emitted 
by the resulting thermal expansion. Acoustic signals are then used to construct 
the deposited thermal energy map. The latter depends on the constitutive optical 
parameters in a nontrivial manner. In this paper, we develop and use an inverse 
transport theory with internal measurements to extract information on the optical 
coefHcients from knowledge of the deposited thermal energy map. We consider 
the multi-measurement setting in which many electromagnetic radiation patterns 
are used to probe the domain of interest. By developing an expansion of the 
measurement operator into singular components, we show that the spatial variations 
of the intrinsic attenuation and the scattering coefficients may be reconstructed. We 
also reconstruct coefficients describing anisotropic scattering of photons, such as the 
anisotropy coefficient g{x) in a Henyey-Greenstein phase function model. Finally, we 
derive stability estimates for the reconstructions. 



Keywords. Photoacoustics, Optoacoustics, transport equation, inverse problems, 
stability estimates, internal measurements. 

1. Introduction 

Photoacoustic imaging is a recent medical imaging technique combining the large 
contrast between healthy and unhealthy tissues of their optical parameters with the 
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high spatial resolution of acoustic (ultrasonic) waves. Electromagnetic radiation, sent 
through a domain of interest, generates some heating and a resulting thermal expansion 
of the underlying tissues. The mechanical displacement of the tissues generates acoustic 
waves, which then propagate through the medium and are recorded by an array of 
detectors (ultrasound transducers). The photoacoustic effect is now being actively 
investigated for its promising applications in medical imaging. We refer the reader 
to e.g. [TOl [131 [IS] for recent reviews on the practical and theoretical aspects of the 
method. 

In an idealized setting revisited below, the electromagnetic source is a very short 
pulse that propagates through the domain at a scale faster than that of the acoustic 
waves. The measured acoustic signals may then be seen as being emitted by unknown 
initial conditions. A first step in the inversion thus consists in reconstructing this initial 
condition by solving an inverse source problem for a wave equation. This inversion is 
relatively simple when the sound speed is constant and full measurements are available. 
It becomes much more challenging when only partial measurements are available and 
the sound speed is not constant; see e.g. [H [71 [TOl [131 [El- 

A second step consists of analyzing the initial condition reconstructed in the first 
step and extracting information about the optical coefficients of the domain of interest. 
The second step is much less studied. The energy deposited by the radiation is given by 
the product of aa{x), the attenuation in the tissue and of I{x), the radiation intensity. 
The question is therefore what information on the medium may be extracted from (Tf,/. 
The product can be plotted as a proxy for aa when I is more or less uniform. This, 
however, generates image distortions as has been reported e.g. in [11]. The extraction 
of e.g. (Ta from aal remains an essentially unsolved problem. 

Two different regimes of radiation propagation should then be considered. 
In thermoacoustic tomography (TAT), low frequency (radio- frequency) waves with 
wavelengths much larger than the domain of interest, are being used. We do not consider 
this modality here. Rather, we assume that high frequency radiation is generated in the 
near infra red (NIR) spectrum. NIR photons have the advantage that they propagate 
over fairly large distances before being absorbed. Moreover, their absorption properties 
have a very large contrast between healthy and cancerous tissues. In this regime, I{x) 
may be interpreted as a spatial density of photons propagating in the domain of interest. 
The density of photons is then modeled by a transport equation that accounts for photon 
propagation, absorption, and scattering; see (12. lip below. 

This paper concerns the reconstruction of the optical parameters in a steady- 
state transport equation from knowledge of H{x) := aa{x)I{x). The derivation of the 
transport equation given in (12. lip below is addressed in section [2l We are concerned 
here with the setting of measurements of H{x) for different radiation patterns. Our 
most general measurement operator A is then the operator which to arbitrary radiation 
patterns at the domain's boundary maps the deposited energy H{x). 

We analyze the reconstruction of the absorption and scattering properties of the 
photons from knowledge of A. The main tool used in the analysis is the decomposition 
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of A into singular components, in a spirit very similar to what was done e.g. in 
[H El E] (see also [2]) in the presence of boundary measurements rather than internal 
measurements. The most singular component is related to the ballistic photons. We 
show that the analysis of that component allows us to reconstruct the attenuation 
coefficient aa and the spatial component as of the scattering coefficient. The anisotropic 
behavior of scattering is partially determined by the second most singular term in A, 
which accounts for photons having scattered only once in the domain. Although the full 
phase function of the scattering coefficient cannot be reconstructed with the techniques 
described in this paper, we show that the anisotropy coefficient g{x) that appears in the 
classical Henyey-Greenstein phase function is uniquely determined by the measurements 
in spatial dimensions n = 2 and n = 3. Moreover, all the parameters that can be 
reconstructed are obtained with Holder-type stability. We present the stability results 
in detail. 

When scattering is very large, then radiation is best modeled by a diffusion 
equation characterized by two unknown coefficients, the diffusion coefficient D{x) and 
the attenuation coefficient Caix). This regime is briefly mentioned in section \2M 

The rest of the paper is structured as follows. Section [2] is devoted to the derivation 
of the stationary inverse transport problem starting from the transient equation for 
the short electromagnetic pulse. The main uniqueness and stability results are also 
presented in detail in this section. The derivation of the uniqueness and stability results 
is postponed to the technical sections [3] and HI The former section is devoted to the 
decomposition of the albedo operator into singular components. Useful results on the 
transport equation are also recalled. The latter section presents in detail the proofs of 
the stability results given in section O 

2. Derivation and main results 

2.1. Transport and inverse wave problem 

The propagation of radiation is modeled by the following radiative transfer equation 



We assume here that S(t, x, v) is compactly supported in t > and in x outside of a 
bounded domain of interest X we wish to probe. The domain X is assumed to be an 
open subset of M" with boundary. For the sake of simplicity, X is also assumed to be 
convex although all uniqueness and stability results in this paper remain valid when X 
is not convex. Here, c is light speed, S*^"^ is the unit sphere in M" and T is the transport 
operator defined as 



where a{x,v) is the total attenuation coefficient and k{x,v',v) is the scattering 
coefficient. Both coefficients are assumed to be non-negative and bounded (by a constant 




(2.1) 




(2.2) 
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M < oo) throughout the paper. We define 

as{x,v) = / k{x,v,v')dv\ aa{x,v) = a{x,v) — as{x,v). (2.3) 

We also refer to scattering coefficient and define aa as the intrinsic attenuation 

coefficient. We assume that for (almost) all (x, t;) G M" x we have aa{x, v) > ao > 
0. We also assume that k{x, v, v') = for (almost) all x ^ X. 

The optical coefficients k{x, v', v) and a{x, v) are the unknown coefficients inside 
X that we would like to reconstruct by probing the domain X by radiation modeled 
by S{t,x,v). In photo-acoustics, the emitted radiation generates some heating inside 
the domain X. Heating then causes some dilation, which mechanically induces acoustic 
waves. Such acoustic waves are measured at the boundary of the domain X. After 
time reversion, the latter measurements allow us to infer the intensity of the source of 
heating. This gives us internal measurements of the solution u of the transport equation 
(12. ip . The objective of this paper is to understand which parts of the optical parameters 
may be reconstructed from such information and with which stability. 

Before doing so, we need an accurate description of the propagation of the acoustic 
waves generated by the radiative heating. The proper model for the acoustic pressure 
is given by the following wave equation (see e.g. [9]) 

np{t,x)=P^H{t,x), (2.4) 
where □ is the d'Alembertian defined as 

□P=^|f-AP. (2-5) 

with Cs the sound speed, where /3 is a coupling coefficient assumed to be constant and 
known, and where H{t, x) is the thermal energy deposited by the radiation given by 

H(t,x)= / aa{x,v')u{t, x,v')dv' . (2.6) 

Not surprisingly, the amount of heating generated by radiation is proportional to the 
amount of radiation u and to the rate of (intrinsic) absorption cTq. 

As it stands, the problem of the reconstruction of the source term H[t, x) inside 
X from measurements of p{t, x) on the boundary dX is ill-posed, because H is 
{n + 1)— dimensional whereas information on (t, x) G M_|_ x dX is n— dimensional. What 
allows us to simplify the inverse acoustic problem is the difference of time scales between 
the sound speed Cg and the light speed c. 

To simplify the analysis, we assume that Cg is constant and rescale time so that 
Cg = 1. Then c in (12.11) is replaced by ^, which in water is of order 2.3 10^/1.5 10^ ~ 
1.5 10^ := ^ ^ 1. The transport scale is therefore considerably faster than the acoustic 
scale. As a consequence, when the radiation source term S{t, x, v) is supported on a 
scale much faster than the acoustic scale, then u is also supported on a scale much faster 
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than the acoustic scale and as a result H{t, x) can be approximated by a source term 
supported at t = 0. 

More specifically, let us assume that the source of radiation is defined at the scale 
e = ^ so that S is replaced by 

Se{t,x,v) = ^p{^)So{x,v), (2.7) 

where p > is a function compactly supported in t G (0, oo) such that Jj^^ p{t)dt = 1. 
The transport solution then solves 

1 d 

-—Ue{t, X, V) + TUe{t, X, v) = Se{t, X, v) . 

With Cs = 1, we verify that is given by 

Ue{t,x,v) = ^u[^,x,v), (2.8) 

where u solves (12.11) with S{t,x,v) = p{t)So{x,v). Because cr^ > (Jq > 0, we verify that 
u decays exponentially in time. This shows that Us lives at the time scale e so that 
H^{t, x) = /g„-i Caix, v)ue{t, X, v)dv is also primarily supported in the vicinity of t = 0. 

Let us formally derive the equation satisfied by Pe when e — > 0. Let (p{t,x) be a 
test function and define (■, ■) as the standard inner product on R x R". Then we find 
that 

{Dp,, ^) = (3{^, = -P{He. ^) = ^(^t))), 

where H{t, x) = /g„-i o"a(x, v)u{t, x, v)dv with u defined in (12.81) and where we have used 
the change of variables t — > et. The latter term is therefore equal to 

-(3{H, ^(0)) = -/3^(0) [ aa{x, v)u{t, x, v)dtdv, 



dt^ " '■^ dt 



up to a small term for sufficiently smooth. We thus find that Pe converges weakly as 
e ^ to the solution p of the following wave equation 

□p = t > 0, X e 

p{{i,x) = Hq{x) := I aa{x,v)u{t, x,v)dtdv, x G C2 g) 

^(0,x) = xgM". 
dt^ ^ 

The inverse problem for the wave equation is now well-posed. The objective is to 
reconstruct Ho{x) for x G X from measurements of p{t,x) for t > and x G dX. Such 
an inverse problem has been extensively studied in the literature. We refer the reader 
to e.g. [13] for a recent inversion with (known) variable sound speed. 

In this paper, we assume that Hq{x) has been reconstructed accurately as a 
functional of the radiation source S'o(x, w). Our objective is to understand which parts 
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of the optical parameters cr(x, v) and k{x, v', v) can be reconstructed from knowledge of 
Hq{x) for a given set of radiations 5*0 (x, v). Note that we will allow ourselves to generate 
many such Sq{x^v) and thus consider a multi- measurement setting. 

The time average u{x,v) := fu(t,x,v)dt satisfies a closed-form steady state 
transport equation given by 

Tu = So{x, v), (x, v)eWx (2.10) 

as can be seen by averaging (12. ip in time since J p{t)dt = 1. 

2.2. Inverse transport with internal measurements 

Since So{x,v) is assumed to be supported outside of the domain X and scattering 
k{x, v', v) = outside of X, the above transport equation may be replaced by a boundary 
value problem of the form 



V ■ + a{x,v)u — / k{x,v' ,v)u{x,v')dv' = 0, {x,v) & X x Ei"- ^ 

i§n-i (2-11) 

u{x,v) = 4>{x,v) (x, f) G r_, 

where the sets of outgoing and incoming boundary radiations are given by 

r± = {{x, v)edX X §""\ ±v ■ u{x) > 0}, (2.12) 

where z/(x) is the outward normal to X at x G dX and 0(x, v) is an appropriate set of 
incoming radiation conditions obtained by solving v ■ V xU + o'a{x,v)u = Sq outside of 
X assuming that aa{x,v) is known outside of X. 

We are now ready to state the inverse transport problem of interest in this 
paper. It is well known that ( 12. lip admits a unique solution in L^{X x §"~^) when 
(p{x,v) G L^(T^,d^), where d^ = \v ■ h'{x)\dp{x)dv with dp the surface measure on dX. 
We thus define the albedo operator as 

A: L\T^,dO - L\X) 

f (2.13) 

0(x, f) I— > A(f){x) = H{x) := / aa{x,v)u{x,v)dv. 



n-l 



The inverse transport problem with angularly averaged internal measurements thus 
consists of understanding what can be reconstructed from the optical parameters a{x, v) 
and k{x, v', v) from complete or partial knowledge of the albedo operator A. We also 
wish to understand the stability of such reconstructions. 



2.3. Albedo operator and decomposition 

The inverse transport problem and its stability properties are solved by looking at a 
decomposition of the albedo operator into singular components. Let a{x,x',v') be the 
Schwartz kernel of the albedo operator A, i.e., the distribution such that 

A(t){x) = j a{x,x',v')(l){x',v')dp{x')dv'. (2.14) 
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The kernel a{x,x',v') corresponds to measurements of H{x) at x G X for a radiation 
condition concentrated at x' G dX and propagating with direction v' G §"~^. Such a 
kernel can thus be obtained as a limit of physical experiments with sources concentrated 
in the vicinity of {x',v') and detectors concentrated in the vicinity of x. 

The kernel a{x, x', v') accounts for radiation propagation inside X, including all 
orders of scattering of the radiation with the underlying structure. It turns out 
that we can extract from a{x,x',v') singular components that are not affected by 
multiple scattering. Such singular components provide useful information on the optical 
coefficients. Let us define the ballistic part of transport as the solution of 

V ■VxUo + a{x,v)uo = 0, in X x §"~\ uq = (f), on r_. (2.15) 

Then for m > 1, we define iteratively 

V ■ VxUm + cr{x,v)um = J k{x, v' , v)um-i{x, v')dv' , in X X Um = 0, on r_. 

(2.16) 

This allows us to decompose the albedo operator as 

A = Ao + A^ + g^, (2.17) 

where Ak for k = 0,1 are defined as A in ( 12.13P with u replaced by Uk and where Q2 is 
defined a.s A — Aq — Ai. Thus, Aq is the contribution in A of particles that have not 
scattered at all with the underlying structure while Ai is the contribution of particles 
that have scattered exactly once. 

Let Ok for /c = 0, 1 be the Schwartz kernel of A^ and r2 the Schwartz kernel of 
Q2 using the same convention as in (12.141) . We define t±{x,v) for x G X and v G S""^ 
as T±{x,v) = inf{s G ]R+|x ± sv ^ X}. Thus, t±{x,v) indicates the time of escape 
from X of a particle at x moving in direction ±v. On dX, we also define 5{x}iy) as 
the distribution such that J^^ 6{x}{y)4>i.y)dlJ>{y) = 4>{x) for any continuous function (f) 
on dX. Finally, we define the following terms that quantify attenuation. We define the 
function E{xo,xi) on X x dX as 

E(xo, xi) = exp i — / cr(xo — S-j r, -j r)(is). (2.18) 

V \xo-xi\ |xo-xi|' / 

We still denote by E the function defined above for xi in X. Then by induction on m, 
we define 

E{xi, Xm) = E{xi, . . . , Xm-l)E{Xm-l, X^). (2.19) 

The latter term measures the attenuation along the broken path [xi, . . . , Xm]- 
Then we have the following result. 
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Theorem 2.1 Let ao, a\, and T2 be the Schwartz kernels defined as above. Then we 
have: 

ao{x,x',v') = aa{x,v')exp a{x - sv' ,v')dsj6{r,-r_{x,v')v'}{x') 



ai{x,x',v') = \iy{x') ■ v' 







aaix,v)-, ; — — -kiX + tf ,f ,:-:,'-t'v' dt . 



\x — x'—t'v' 



(2.20) 

Moreover, we have the bound 



\x — x' — {{x — 


x') ■ v')v' 


|"-3r2(x,x',i;') 




z/(x') • v'\ 





r2(x, x', (,Loo^XxT_) when n = 2 

(2.21) 

e L°°(X X r_) when n>3. 



Ivix') ■ v' 



This theorem wiU be proved in section [3l The resuhs show that is more singular 
than «! and It turns out that cti is also more singular than r2 in the sense that it 
is asymptotically much larger than r2 in the vicinity of the support of the ballistic part 
Oq. That this is the case is the object of the following result (see also Lemma [3.21 (13.91) 
and (13.111) for "m = 1"). For any topological space Y, we denote by Cb{Y) the set of 
the bounded continuous functions from Y to M. 

Theorem 2.2 Let us assume that a e Cb{X x §"~i) and that k e Cb{X x §"-1 x §"-i). 
For (x',t/) G r_, let x = x' + tpf' for some G (0, r+(x',f')). Let v'-^ G §"~^ be such 
that v' ■ f'-*" = 0. Then we have the following asymptotic expansion: 

aiix + ev'^,x' ,v') /, 1\/ / , ,n / , ,n\ 1\ 

— = \yi-)[x[x,v ,v) + x[x,v ,-v)) +0 In- 

E[x,x')\v[x') ■ v'\ \ eJ ^ ' \ e/ (2 22) 



sm 



n— 3 



q;i(x + , X , f 

E(X.X')\U{X')-V'\ £"-2/ — ^ , ^v"//-" . -y^. 



ex{.x,v\v{e))de + o 



for n = 2 and n > 3, respectively, where we have defined the functions x{x,v',v) = 
(Ta(x, v)k{x, v', v) and v{9) = cos 6v' + sin 6v'^ . 

Theorem 12.21 will be proved in section [3J 

We thus observe that ai blows up a priori faster than r2 as e ^ 0, i.e., as 
the observation point x becomes closer to the segment where the ballistic term ao 
is supported. This singularity allows us to obtain information on the optical coefficients 
that is not contained in the ballistic part ao. Moreover, because of the singular behavior 
of «!, such information can be reconstructed in a stable manner. 



2.4- Stability estimates 

As mentioned above, the singular behaviors of and ai allow us to extract them 
from the full measurements a. Moreover, such an extraction can be carried out in a 
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stable fashion, in the sense that small errors in the measurement of the albedo operator 
translates into small errors in the extraction of the terms characterizing and ai. 

More precisely, let A be the albedo operator corresponding to optical parameters 
(cr, k) and A the operator corresponding to the optical parameters (a, k). From now on, a 
term superimposed with the ~ sign means a term calculated using the optical parameters 
(o", k) instead of (a, k). For instance E{x, y) is the equivalent of E{x, y) defined in (12.181) 
with (a, k) replaced by (a, k). 

We first derive the stability of useful functionals of the optical parameters in terms 
of errors made on the measurements. Let us assume that A is the "real" albedo operator 
and that A is the "measured" operator. We want to obtain error estimates on the useful 
functionals of the optical parameters in terms of appropriate metrics for A — A. We 
obtain the following two results. The first result pertains to the stability of the ballistic 
term in the albedo operator: 

Theorem 2.3 Let A and A be two albedo operators and {x',v') G r_. Then we obtain 
that 



r+{x',v') 



< 







A\ 



dt 



(2.23) 



Theorem 12.31 is proved in section 14. 1[ 

The stability result obtained from the single scattering component is based on a 
singular behavior obtained in the vicinity of the ballistic component. Such a behavior 
cannot be captured by the norm used above. Instead, we define Fi = a — ao as the 
Schwartz kernel of the albedo operator where the ballistic part has been removed, i.e., 
for measurements that are performed away from the support of the ballistic part. Our 
stability results are obtained in terms of errors on Fi rather than on A. We can then 
show the following stability result. 



Theorem 2.4 Let us assume that {a, a) G Cb{X x §" ^)^ 
^ gn-i)2^ (^^^^/) ^ X xdX and define v' - 



V ■ V 



< 



0. In dimension n = 2, we have 



\x—x'\ 



and that {k,k) G Cb{X x 
Let v'^ G S""^ such that 



E{x, x'){x{x, v', v') + x{x, v', -v')) - E{x, x'){x{x, v', v') + x{x, v' , -v')) 
{Ti-Ti){x,x',v') 

lulx') ■ V'\W2{X,X' ,V') L°°(Xxr_)' 



(2.24) 



where W2{x, x' , v') = 1 + In ^1^—2 — '^^ ) ■ When n > 3, we h 



sin^-%e) E{x, x')x{x, v', v{e)) - E{x, x')x{x, v' , v{e)) dd 



fFi -VMx.x'.v'^ 



ave 



(2.25) 
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where Wn{x,x' ,v') = \x — x' — ((x — x') ■ v')v'\'^ and where we use the same notation 
as in Theorem \2.B. 

Theorem 12.41 is proved in section I4.2[ 

Such results do not grant uniqueness of the reconstruction of the optical parameters 
in the most general setting. However, they do provide stable, unique, reconstructions 
in several settings of interest. 



2.5. Scattering- free setting 

Let us first assume that /c = so that a = ctq. Then knowledge of the albedo operator 
uniquely determines aa{x,v) for all x G X and v G S*^"^. 
Indeed, we deduce from Theorem 12.31 that 

is uniquely determined and hence e~ ca{x'+sv' ,v')ds gjj^pg ^y^q latter equals 1 when t = 0. 
Taking the derivative of the negative of the logarithm of the latter expression gives us 
(j{x' + tv',v') for all {x',v') G r_ and t > and hence a{x,v) for all {x,v) G X x S"~^. 
Moreover, we have the following stability result. 

Theorem 2.5 Recalling that aa{x,v) is bounded from above and below by positive 
constants, we find that when k = 0, 

Ika - 5-a||Loo(gn-i.£^i(jf)) < C\\A - A\\c(L^(^r.,d^y,L^{x))- (2.26) 
Here, C is a constant that depends on the uniform bound M. 

The above theorem is proved in section 14.11 Note that the above result is local in x' 
and v'. In other words, o"a(x, v') is uniquely determined by {a{y, x — r_(x, v')v', v')\y = 
x + tv' for — r_(x, v') < t < r+(x, v')}, that is by the experiment that consists of sending 
a beam of radiation in direction v' passing through the point x (at least asymptotically 
since such a transport solution is not an element in L^{X x 



2. 6. Reconstruction of the spatial optical parameters 

We now assume that k ^ Q. Then the ballistic component of A and the estimate 
in Theorem 12.31 allow us to uniquely reconstruct both aa{x,v) and a{x,v) under the 
assumption that 

aa{x,v) = aa{x,-v), a{x,v) = a{x,-v). (2.27) 

We recall that o"a(x, v) is bounded from below by ctq > 0. For technical reasons, we 
also assume that a{x, v) is known in the ^o^vicinity of dX, i.e., for all (x, v) G M" x S"~^ 
such that dist(x,9X) < 6o for some Sq > 0. Such a hypothesis is not very restrictive 
from a practical viewpoint. 
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We denote by W ^'^{X) the Banach space of the continuous hnear functionals on 
the Banach space Wq'°^{X) := {(/> E L~(X) | suppcp C X, V0 G L°°(X,C'^)} (where 
the gradient V is understood in the distributional sense). 

Under the above assumptions, we have the following result. 

Theorem 2.6 The coefficients aa{x,v) and cr(x, f) are uniquely determined by the 
albedo operator A. Moreover, we have the following stability estimate 

Ik - cr||Loo(§n.i.vi/-.i,i(x)) ^ ^2 28) 

+ ||0"a — Cra||L°°(§"-i;Li(^)) ^ C\\A — A\\c[Ll(^^^,d^)■,L^(X))■ 
The constant C depends on the parameter Sq as well as the uniform bounds ctq and M. 

This theorem is proved in section |4?T1 The above result is also local in x' and v'. In other 
words, aa{x,v') and a{x,v') are uniquely determined by {a{y,x — t^{x,v')v' ,v')\y = 
X + tv' for — r_(x, f') < t < r+(x, f')}, i.e., by the experiment that consists of sending 
a beam of radiation in direction v' passing through the point x. 

Stability estimates may be obtained in stronger norms for a provided that a priori 
regularity assumptions be imposed. We show the 

Corollary 2.7 Let us assume that a and b are bounded in L°°{E>"'~^,W'^'''^{X)) by Cq 
forp > 1 and r > —1. Then for all —1 < s < r , we have 

\\a - cr\\La.(sn~i.w='P(x)) < C\\A - ^||£('ii''(r_,d5);Li(x))' (2.29) 

where the constant C depends on Co and on the uniform bounds ctq and M. 
The corollary is proved in section 14. 1[ 

2.7. Application to Henyey-Greenstein kernels 

Let us assume that cr G C6(X x §""!) and that A; G Cb(X x x §"^i) and let us 

assume again that a{x) is known in the 5o~vicinity of dX, i.e., for all x G M" such that 
dist(x, dX) < 6o for some 6o > . 

The stability estimate fl2.23p allows one to uniquely reconstruct cXa and a under 
the symmetry hypothesis (12.271) . which is quite general physically. Indeed, even when 
attenuation is anisotropic, there is no reason to observe different attenuations in 
direction v and direction —v. The stability estimate in Theorem 12.41 provides additional 
information on the optical coefficients, but not enough to fully reconstruct the scattering 
kernel k{x, v',v). 

In dimension n = 2, we gain information only on k{x,v',v') + k{x,v' , —v'). In 
dimension n > 3, we garner information about JJ^ ^(^^i "^'j cos6't'' + sin6v'^)d6 for all 
f'-*- orthogonal to v'. The integration in 6 means that one dimension of information is 
lost in the measurements. Thus, 3n — 3 dimensions of information are available on the 
(3n — 2)— dimensional object k{x,v',v). 
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Let us consider the case of an isotropic absorption 0"^ = <^a{x) and isotropic 
scattering in the sense that k{x, v', v) = k{x, v' -v). In such a setting, k{x, v' ■ v) becomes 
(n + 1)— dimensionaL Yet, available data in dimension n >2 give us information on 

ag{x) := k{x, -1) + k{x, 1), when n = 2, (2.30) 

ag{x) := / A;(a;,cose)sin"~^ede, when 72 > 3. (2.31) 
Jo 

This is different information from the normalization in (12. 3p 



a.Ax 



[ k{x,v' ■v)dv' =\W^\ f k{x, COS 9) sin""-^ Ode. 

J§"-i Jo 



As a consequence, if k{x,cos6) is of the form cis (x)/(x, cos 6^), where f{x,cos9) is 
parameterized by one function g{x), then we have a chance of reconstructing g{x) from 
knowledge of o"g(x) and as{x) provided as{x) > (where = c — aa). 

This occurs for the classical Henyey-Greenstein (HG) phase function in dimensions 
n = 2 and = 3, where 

k{x, A) := ^^(^) 2,(i^'^^^'/'^2,(x)A) ' " = 

1 - g'^ix) 

k{x, A) := c^.(a;)-— — — when n > 3, (2.33) 

47r(l + 5((x)2 - 2g{x)\)2 

where g G Cb{X) and < 5'(x) < 1 for a.e. x E X. Note that 

ag{x) = as{x)h{g{x)), (2.34) 



for X E X where the function /i : [0, 1) ^ M is given by 

1 + ^^ 
7r(l — K 



h{K) := — --, when n = 2, (2.35) 

^ / 1 --2 \ 



hU) := / ^ci^, when n > 3. (2.36) 

Jo 47r(l + ^2 - 2fi;cos(6'))2 

Theorem 2.8 In the HG phase function in dimension n = 2,3, the parameter g{x) is 
uniquely determined by the data provided crs(x) > for a.e. x E X . 

Theorem ES] follows from Theorems EJl El and from flOOD . flCTl) . flOD and 
the following Lemma. 

Lemma 2.9 The function h is strictly increasing on [0,1), h{0) = 0, h{0) > 0, 
limg^i-(l — g)h{g) = c{n), where h{g) = ^^{g) and c(2) = vr"^ and c(3) = {2it)~^. 

Considering (12.351) . Lemma [2.91 in dimension n = 2 is trivial. The proof of Lemma [2.91 
in dimension n = 3 is given in section 14.21 

Moreover we have the following stability estimates. 
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Theorem 2.10 In dimension n>2, we have 

(ri-fi)(x,a;',i;') 



\l>{x') ■ v'\'Wn{x, X', V') 



+ 11^-^1 



C(L\T-,d£,);L^X))]: 



(2.37) 

where Wn{x,x' ,v') is defined in Theorem \2.4\ and where the constant C depends on 6o, 
(To and M . 

In addition, for the HG phase function in dimension = 2, 3, we have the following 
stability estimate: assume mhi{as[x),as{x)) > cr^.o > and {h{g{x)), h{g{x))) G 
W^'°°{Xf, then 



\\h{g{x))-hCg{x))h^^x)<C 



\A-A\ 



{T,-f^){x,x',v') 

\c{L^r.,diy,LHX))J, 
(2.38) 

where the constantC depends on6o, ctq, asfl, M andYmYi{\\h{g)\\-^i,oo(^x)^ \\h{.g)\\w^-'^{x)) ■ 

Theorem 12.101 is proved in section 14.21 

The sensitivity of the reconstruction of g{x) degrades as g converges to in the 
sense that C{Go) — oo as Go — ^ 0, Go = max(||5f||oo, ||^||oo), where G is the constant 
that appears on the right-hand side of f l2.38p when we replace h{g{x)) — h{g{x)) by 
g — g on the left-hand side of (12.381) . On the other hand, C{go) when go — > 1, 
go = min(||5f||oo, ||^||oo), so that reconstructions of g{x) are very accurate for g{x) close 
to 1, i.e., in the case of very anisotropic media. 

More precisely, using the properties of the function h, one can replace the left- 
hand side of (12.381) by \\g^ - g'^Wmx) (resp. H^f - resp. \\^ - j^\\l^(x)) 
provided that max(||(7||oo, ||^||oo) < Go < 1 (resp. < go < min(||5f||oo, ||^||oo) and 
max(||5f||oo, ll^lloo) < Go < 1, resp. min(||5f||oo, ||^||oo) > fl'o > 0) for some constant Go 
(resp. ((y'o,Go), resp. go) and the constant G on the right hand side of (I2.38P then 
depends also on Go (resp. {go, Go), resp. go). 



2.8. Reconstructions in the diffusive regime. 

When scattering is large so that the mean free path ^(x) is small and intrinsic 
attenuation (Ja{x) is small, then radiation inside the domain X is best modeled by 
a diffusion equation 

-V-D{x)VI{x)+aa{x)I{x)=Q xeX 39) 
I{x) = (j){x) X G (9X, 

where /(x) = fg„_j^u{x,v)dv is the spatial density of photons and D{x) is the 
diffusion coefficient. We refer the reader to e.g. [2[ [8] for references on the diffusion 
approximation. When scattering is e.g. isotropic, i.e., when k{x,6',6) = k{x), then we 
find that D{x) = ^J'^^^^ , where ag is introduced in (12.31) and n is the spatial dimension. 

When D{x) is known, then the reconstruction of o"a(x) may be easily obtained 
by using only one measurement. Indeed, the measurement H{x) = aa{x)I{x) so that 
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J(x) may be obtained by solving (12.391) . Once I{x) is known, it will be positive in X 
provided that </'(x) is non trivial and non- negative. Then o"a(x) is obtained by dividing 
H by /. When (t){x) is bounded from below by a positive constant, then we see that the 
reconstruction of (Tq is unique and clearly stable. 

When (D(a;), (Ta(x)) are both unknown, then multiple (at least two) measurements 
are necessary. This problem will be analyzed elsewhere. 

This concludes the section on the derivation and the display of the main results. 
The mathematical proofs are presented in the following two sections. 



3. Transport equation and estimates 

In this section, we prove several results on the decomposition of the albedo operator 
(Lemmas 13. H 13.21 and 13.31) and prove Theorems 12.11 and 12.21 

We first recall the well-posedness of the boundary value problem (12. lip (here 
a G L°°(X X S"-i), keL°°{Xx §"-1 X §"-i) and > (Tq > 0) and give a decomposition 
of the albedo operator. The boundary value problem (12.111) is equivalent to the integral 
equation 

{I-K)u = J(t) (3.1) 

for u G L^{X X §"~^) and G L^(r^,dC,), where K is the bounded operator in 
L\X X §"-!) defined by 

/■T-{x,v) r 

Ku = I E{x,x — tv) / k{x — tv,v' ,v)u{x — tv,v')dv'dt, (3.2) 

Jo JS'i-i 

for a.e. {x,v) G X x S""^ and for u G L^{X x §"~^), and J is the bounded operator 
from L\T_,d^) to L\X x §"-i) defined by 

J(f){x,v) =e-^o"*'''"'"(^-^"^'^)'^P0(x-r_(x,w)i;,w), (3.3) 

for a.e. {x,v) G X x §""1. 

Since = a — cTs > 0, it turns out that J — is invertible in L^{X x S'"^^) [31 
Lemma 2.4], so that the albedo operator A : L^{r^,d^) L^i^) is well-defined by 
(I2.13P and so that the solution u of (12. lip with boundary condition (p E (r__ , d^) 
satisfies 



n 

U = 

m.=0 

It follows that 



^i^"'J(/> + i^"+^(/-7r)-V(/). (3.4) 



A(p = / j 'y^am{x,x' ,v') j (f){x' ,v')dfj,{x')dv' 
■^r- Vm=o / 



'jn+i{x,y,w){{I - K) ^J(f)){y,w)dydw, (3.5) 
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where 7^, m > 0, is the distributional kernel of : L^{X x S"~^) L^{X) defined 
by 

K^u{x)= / aaix,v)K"'u{x,v)dv, (3.6) 



for a.e. x G X and u G L^{X x S" and where am, m > 0, is the distributional kernel 
of if"* J : Li(r_, dO L\X). 

We give the explicit expression of the distributional kernels a^, m > 2, and 7^, 
m > 3 in Lemma 13.11 and study the boundedness of am in Lemma 13.21 and of 7n+i in 
Lemma 13. 3[ We then prove Theorems 12.11 and 12.21 and Lemmas 13.21 and 13. 3[ The proof 
of Lemma [3.11 is given in the appendix. For w G M", w 0, we set w := j^. 

Lemma 3.1 Form > 2 and a.e. (zq, Zm,Vm) G X x r_, we have 

W{Zm)-Vm\ I I [(^a{Zo,Vo) (3.7) 

X 



Jo 

-^(^-O) • • • ) •2'm— 1) •2'm ~l~ t Vmi Z^m)k(^Zm ~\~ t Vmi Vm: '^rn—l) 



1 7 -Lf'v —7 Jn-lTTm-2| _ |n+l 
^nj=l ^(^^j) I'j) "^i-l)] tii^Zi^^^i, i=o...m-2, C^^l • • • dZm-1. 



2j^l. i— 0...m — 2, 



For m > 3 anc? a.e. {zq, Zm, Vm) G X x X x ^ , we have 

-fm{zo,Zm,Vm) = ^^^^ . (3.8) 

X [aaizo, Vo)UZik{zi, Vi, Wi-i)]„^=,— ^, i=o...m-i dzi... dZm-l- 
Lemma 3.2 For n = 2, 

"^^"'"'"'^ ^GL-(Xxr_), (3.9) 



\iy(x') ■ v'\ m / n +^ ' ^ ^ 

I ^ ' I \ \x — x'\ — (x — x')-v 

##4 e ^^(^ X r-). (3.10) 



For n >3, we have 
\x — x' — {{x — x') ■ v')v'\'^~^~'^am{x, x', v') 



G L°°{X X r_), 1 < m < n - 2, (3.11) 



\v{x') ■ w'l In 



an-i{x,x\v') eL°°(Xxr_), (3.12) 



\x—x'— {{x—x')-v')v'\ 



""I'^'f'';!^ G L-(X X r_). (3.13) 
|z^(x') ■ v'\ 

We do not use (13.91) and (13. lip for m = 1 in order to prove (12.211) . However we will 
use them in the proof of the stability estimates given in Theorem 12. 4[ 
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Lemma 3.3 For n >2, we have 

7„+i G X X X (3.14) 

In addition, we have 

K^'+^il -K)-^J(t){x)= I Tn+i{x,x',v')(l){x',v')dfx{x')dv', (3.15) 



for a.e. x G X and for (j) G L^(r_,(i^), where 

e L°^(x X r_). (3.16) 

\i^[x') ■ v'\ 

Proof of Theorem \2.1i The equality fl2.20p for ao follows from the definition of the 
operator J ([SJl). From ([MD, ([S^D and it follows that 

KJ(p{x) = / aa{x,v) / / E{x, X — tv, X — tv — T^{x — tv,v')v') 

X /c(x — tv, v' , v)(f>{x — tv — T-{x — tv, v')v', v')dv'dtdv, (3-17) 

for a.e. x & X and (p G L^{r^,dC,)- Performing the change of variables z = x — tv 
(dz = f^-^dtdv, t = I z — x\) on the right-hand side of ( 13.17p . we obtain 

KJm= [ ^f^^^^^E{x,z,z-T.{z,v')v') 

JXx§"-i \X — Z\ 

X k{z, v' , X — z)(j){z — r_ {z, v')v', v')dv'dz, (3.18) 

for a.e. x G X and (j) G L^{r^,d^). Performing the change of variables z = x' + tv' 
(x' G dX, t > 0, dz = |z/(x') ■ v'\dtdfi{x')) on the right hand side of fl3.18p . we obtain 

KJ^ix)= I " ^ ^ ~ -'y^ Fix, x' + t'v', x') 

Jt- Jo \x — x' — t'v'l"- ^ 

X k{x' + t'v' ,v' , X — x' — t'v')(f>{x' ,v')dtd^{x' ,v'), (3.19) 

for a.e. x G X and (p £ L^{^~, dC,), which yields (12.201) for ai. 

Now set r2 := J2m=2^-rn + Pn+i when n > 2. Taking account of Lemma 13.21 
(I3:T0|) - (I37[3|) and Lemma [^31 (13361) . we obtain ^M>- □ 

Proof of TheoremlKM We assume that a G Ch{X x §""!) and k G Cb{X x §""1 x §""!). 
Let (x',t/) G r_ and let t'^ G (0, r+(x', f')) and let w'"*" be such that v' ■ t;'"*" = 0. Set 
x = x' + t'^v'. From it follows that 



ai(x + 5f x', f ') ^ X, X , / X 

^ ' ' - = (L+ + L_)(5), where (3.20) 



|z/(x') ■ f' 



L+{e) := I ^ ' ^ Zy ^ aa(x + ev'^, Vt',e)k{x' + t'v', v', Vt',s)dt', (3.21) 

Jo {{t' - t'o)^ + e^) 2 



L_{e) := / ^— -i -— ^a„(a; + ev'^ , v,^,)k{x' + tV, v\ v,,e)dt\ 
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{x + ev'~^, x' + t'v', 

(3.22) 

for e e (0, r+(x, v'^)), where w*/,^ = ^i!^v' + w'^ for t' e M. 

We prove fl2.22l) for n = 2. Consider the function arcsinh : R ^ M defined by 
arcsinh(|/) := ln(y + ^/T+~y^), for ?/ G M. Then performing the change of variables 

rj : = '^''^'^^'^^ s — I ^]-^g right hand side of (13.211) and performing the change of variables 

arcsinh( — ) 

■— s ) ^]-^g right hand side of (I3.22p . we obtain 



arcsmn( — ^ ) 



L±(e) := arcsinh(y)L'i(e), (3.23) 

for e G (0, r+(x, w'"*")), where 

s+:=t'^, s.:=T+{x',v')-t'„ (3.24) 

L'i(e) := / E{x + ev'^,x' + t'^{r],e)v',x') 
Jo 

aa{x + ev'-^, v±{ri, e))k{x' + t'^{ri, e)v', v\ v±{r], e))dri (3.25) 

t'^{r],e) := t'^ =]= esinh (77 arcsinh (—)) , (3.26) 



e 



± sinh (?7arcsinh(^)) v' 



'sinh (?7arcsinh(^)) + 1 

for ?7 G (0, 1) (we recall that sinh(y) = ^^^^f^, y eR). Note that using the definition of 
sinh and arcsinh, we obtain 



l^sinh (r7arcsinh( — ) j | < , (3.28) 

for ?7 G (0, 1). Therefore using (13.261) we obtain 

4(r/,e) ^ t'o, as £ 0+ (3.29) 

for ?7 G (0, 1) and i = 1, 2. Note also that from (13.271) it follows that 

±t/ + sinh (?7 arcsinh ( — )) ^ v'^ 
Vi{r],e) = ^' ^-^^ ±v', (3.30) 

v/l + sinh (?7arcsinh(^)) 

for rj G (0, 1) (we used the limit sinh (?7arcsinh(^)) — > +00 as £ -^^ 0^ which holds for 
any positive real numbers s and 77) . 

Using (13.251) . (I3.29p . (13.301) and continuity and boundedness of a and k and a a, and 
using Lebesgue dominated convergence theorem, we obtain 

L'^is) ^ aa{x,±v')E{x,x')k{x' + t'Qv',v',±v') as £ ^ 0+. (3.31) 
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Finally note that 

arcsinh(^) = In (^) + o(^^) , as e ^ 0+. (3.32) 

Combining ([320]), ([323]), ^M) and ^M), we obtain (12:221) for n = 2. 

We prove (I2.22p for n > 3. Performing the change of variables rj = on the 
right-hand side of (I3.2ip . we obtain 



(3.33) 



for £ G (0, r+(x, f '"*")), where ^3(5, r^) = tg + £77 for 77 G M. Note that 

t^j(e,r/) ^ as e 0+. (3.34) 

Therefore, using the Lebesgue dominated convergence theorem and continuity and 
boundedness of (a, k, aa), we obtain 

ro aa{x,v)k{x,v',v) 
L^{e) = e'--Eix, x') / ^ tT'^ dv + o{e'-)^ (3.35) 



as e ^ 0"*". Similarly performing the change of variables rj = on the right hand 
side of f l3.22p . and using Lebesgue dominated convergence theorem and continuity and 
boundedness of (a, k, aa), we obtain 

L.{e)=e'-'^E{x,x') / , ^'T'^+^ dn + o(e'--). (3.36) 







as e —>■ O"*". Note that performing the change of variables cos(^) = ^y^^ ^ ^ ^ (O;^); 
iv = -W§;'dr] = -^dO) we have 



sm(6») ' ' sm{e) 

aa{x,v)k{x,v\v) 



\v=--nv'+v'^ J _ I -^/mn-S 



dr]= / sm{9)''-^aa{x,v{9))k{x,v',v{9))d9. (3.37) 







Adding fl335]) and fl336D and using fl337l) and ([SD, we obtain (12:221) for n > 3. □ 

Proof of Lemma \3.Si First note that from f l2.20p and (13. 7p . it follows that 

am{x,x\v') < \\aa\\oo\\k\\Z\u{x') ■ v'\Im,nix,x' ,v'), (3.38) 
for a.e. (x, x', v') ^ X x r_ and for m G N, m > 1, where 



■l-+ix',v') 1.1 
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for (^o, x', v') & X X r_ and m > 1. 

We prove ([33]) and fl330D . Let n = 2. Let {x,x',v') G X x r_ be such that 
X 7^ x' + Xv' for any A G M. Set {w)± := w — {w ■ v')v' for any w G M". Using f l3.39p and 
using the equahty \x — x' — i!v'\^ = (t' — (x — x') ■ v'Y + — x')_lP, we obtain 

/tj^{x' ,v') — (x—x')-v' 
^r-rdt, (3.41) 



X - 


-x'- 


- T+(x', — 


(x 


-x' 


-r+(x',^;')^') 


■v' 






X — x' — 


(x 


-x') 


1 ■ v' 





In ^ L J' ' ' ' h (3-42) 



where we used that ~ In (x + vT+~x^) . Estimate f l3.9p follows from (13.381) 

and ([3321). 

Let (x, x', i;') G X X r_ be such that x 7^ x' + \v' for any A G M. Using (13.421) (with 
"x = 2") and ([33QD, we obtain 

( \Z~X' ~T+(x' ,V')V'\~(Z~X' ~T+{x' ,V')V')-V'\ ^ 



l2,n[X,X,V)= / ^ . . -dZ 

Jx \x-z\ 
= J^^Jx, x' + r+(x', v')v', v') - J^^Jx, x', v'), (3.43) 

where 

^/ / /N f ln(|2 — a| — (2; — a) ■ w') , / . 

J ,Jx,a,t;' = / ^ ^ ^rf^, 3.44 

Jj^ \x-z\ 

for a G X. We prove 

sup „(x, a, w')| < cxD. (3.45) 

Then estimate (l310|) follows from ([SJH]), (13:431) and (131451) . 

Let (x, a,v') e X X X X Note that by using (13.441) . we obtain 

IP f 'w ^ n f {\z - a\ + {z - a) ■ v')l 

|/2^„(x, a, )| < C / jdz 

Jx \x — z\{\z — aP — ((z — a) ■ v'Y)» 

< {2D)-^C [ i -dz, (3.46) 

Jx \x — z\{\z — aP — {{z — a) ■ v'y)^ 

where C := sup^,g(Q g£)) rs] ln(r)| < 00 and D denotes the diameter of X. Consider 
unit vector orthogonal to v' and perform the change of variable z = 
a + Ai?;' + \2v'^, then we obtain 

|/^,„(x,a,.')| < (2D)ic / ^^.47) 

ihD,D]2 ((Ai,^ - Ai)2 + (Aa,^; - 2 l^ah 

where Xi^xv' + Aa,!;"?^'"*" = x — a. Then note that 

dXi _^^( /^-Ai,, + ((A2-A2,.)2 + (/^-Ai,,)2)i 



D {{\l,x - Al)2 + (A2,. - X2))- \-D- Ai,, + ((A2 - A2,.)2 + (D + Ai,,)2))2 

(3.48) 
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Combining fl3.47p and (13.481) . we obtain 



\I'{x,a,v')\<C' 



D \X2\H-{D + Xi,.) + x/(A2-A2,.)2 + (D + Ai,,)2)l 



1 X 1 



< (D + Am + ((A2 - A2,.)^ + {D + Ai,.)^)t)i ^^^^ 

-D IA2I4 IA2 - A2,a;|4 



<C""/ " . =C"|A2,.|^ / ' , . (3.49) 



_D IA2I4IA2 - A2,xh ' J^J>^ \^\*V'^ 



where C := (2D)iC2, C" := {12D^)iC^ and C" := {72D^)IC\ Finally note that 

\X2,x\^ I —t;^ ^rfs < Ci(A2,x) + C72(A2,x), where (3.50) 



Ci(A2,x) := |A2,x|W l^r^ls- < £)W IsT^ls - l|"2ds (3.51) 

J|s|<2 "'|s|<2 

C2(A2,x) := |A2,alW kr^ls - < |A2,x|^23 / |s|"^rfs 

J2<lsl<max(2.TT-^) J2<lsl<max('2.TT-^) 



'2<|.|<max{2,p^) J2<|s|<max{2,p^) 

< 2t('|A2,x|max(2, -^))' < 2^L)i (3.52) 

Combining ([S3H])-(IS32]) we obtain (1X151) . 

The statements fl3.lip - fl3.13p follow from fl3.38p and the following statements 
f l333l) -f l335l) 

((x-x') ■?;')^T"^"'"Vn(a;,a:',?;0 G L°°{X x r_), (3.53) 



In 



/„,„(x,x',t;0 eL°^(Xxr_), (3.55) 
for {m,n) gNxN, n>3, l<m<n — 2 and where Jm,„ is defined by f l3.39p and 

dsaoD. 

We prove fl3.53p - fl3.55p . which will complete the proof of Lemma [3.21 We proceed 
by induction on m. We prove fl3.53p for m = 1. Let n > 3. Note that formula (13.411) 
still holds. Note also that 

\w - Xv'\^ = \w±\^ + \wv'- A|2 > 2-Ww±\ + \wv'- A|)^ (3.56) 

for {w,X) e X {-D,D), \w\ < 2D. Therefore 

dX ^ f'^ 2^dX ^ 2^ 
.2z,k-At;'|-^ "y-2D(k±l + |A|)"-^ - (n-2)|^^|«-2- ^ ' > 

Thus fl333D for m = 1 follows from fl337D and ([Ml]). 
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Let m G N, m > 1 be such that fl3.53p - fl3.55p hold for any n > 3. We prove that 
fl3.53p - fl3.55p hold for any n > 3 and for "m"= m + 1. Let {x,x',v') be such that 
x' + \v' for any A G M. From fl339D and flH^iOD it follows that 

W,„(x,x',^')< / Y^''i'_''P dz. (3.58) 



X — z\ 



Assume that m + 1 < n — 1. Then from fl3.58p and fl3.53p for (m, n) it follows that there 
exists a constant C (which does not depend on {x,x',v')) such that 

Im+lnix,x'v') < C [ -. ; —7 7T — ; --, dz. (3.59) 

Performing the change of variables z — x' = z' + Xv' , z' ■ v' = 0, we obtain 



z'\<D 



D 



dX \ dz' 



W,„(x,a:',.') < C /^,^,^^ \x-x'-z'-Xvr-'] W^^' ^'"'^^ 



Combining fl3.60p and (13.570 (with "w = x — x' — 2;"'), we obtain 

lyn+iA^^x' ,v') < -r^^,^{x,x\v'), where (3.61) 

n — 2 ' 



dz' 



z'\<D ^ III 



dz' 



L'eB„-,io,D) Wi^- x')±\ei - z'\--^\z'\-^--^' ^^-^^^ 
and ei = (1,0,..., 0) G M""^ and 5„_i(0,L') denotes the Euclidean ball of W'^ 
of center and radius D. Using spherical coordinates z' = |(x — x')_L|ei + rQ, 
(r, n) G (0, +00) X §"-2 ^ ^ (sin(0), cos(^)e) {{6, 6) G (-f , f ) x S"-^^ Q^^^^j^ 

C+i,„(x,x',t;')<l§""'l /'cos"-3(^) 

^3T3;;r)c^^ (3.63) 

/o (|r+ |(x-xO±|sin(^)|2 + |(x-x')±|2cos(^)2)'^/ 

(by convention := 2). Note that by performing the change of variables "r = 
r + |(x — x')_l| sin(^)" and using the estimate + fe^ > 2~^(a + 6)^, we obtain 

dr 



< 



'0 (|r + |(x - x')±| sin(6')P + |(x - x')±P cos(6')2) — 2 — 

2'^dr ^ /-^^ 2^dr 

-3Z) (kl + cos(^))"-i-'" "Jo (r + |(x - xO±| cos(^))'^-i- 

2^ 

ifm + l<n — 1, 



. (n-m-2)(|(x-xO±|cos(^))"-2— ^ .3 Q4^ 

- ^ n+i /3D + |(x-x')±|cos(^)\ ^ ^ ^^■'^ ^ 

2~ In -—^ . / ' ,„r ' , if m + 1 = n - 1, 



3D + |(x-x')± 




|(x - x')±| 


cos(^) ) 
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for ^ G (— |, |). Assume m + 1 < n — 1. Then combining 03.631] and (13.641] . we obtain 



2^ 


gn-3 




{n — 2 — m)\{x — X 


0± 


n—2—m 



cos"'-\e)de. (3.65) 



Therefore using also ( 13.6ip we obtain that ( 13.53P holds for "m" = m + 1 < n — 1. Assume 
m + 1 = n — 1. Then note that 

AD 



^^. SD+\{x-x')^\cos{e) 



\{x — x')± \ cos{6) 
< ln(4L)) + ln(- ^ 



[X — X )±\ 



\{x — x')± \ cos{6) 
ln(cos(^)), 



for ^ G (-f , f ). Combining (l333|) . fITOl] and (Km^ . we obtain 



(3.66) 



(3.67) 



where Ci := /_\ cos'^"^(^)(ln(4D) - ln(cos(e)))rf^ < oo and C2 := /_\ cos"-3(e)d^. 
Therefore using also (13. 61 1) , we obtain that (I3.54p holds for "m"= m + 1 = n — 1. 

Assume that m + 1 = n. From (13.581) and (13.541) for (n — 1, n) it follows that there 
exists a constant C (which does not depend on {x,x',v')) such that 



In 



1 



|n-l 



'X F ~ ^1 

Performing the change of variables z — x' = z' + \v', z' ■ v' = 0, we obtain 



(3.68) 



z''u'=(J 
z'\<D 



D 



dX 



\x — x' — z' — Af'l" 
Combining (I3.69P and (13.571) (with ''w = x — x' — z'"), we obtain 



dz'. (3.69) 



Jm+i,„(a;,x',w) < ^ _ ^ I'l^{x,x',v') 



(3.70) 



where 



'■i''=0 
\z'\<D 





In 1 






\{x- 


x') 







In 


[\''\) 




ll(^- 









-dz' < 



2'6B„-i(o,D) 2 ||(x — x')_L|ei — 2 



(3.71) 



and C" := sup^,g(o,D) '"^ I liil'^)!; = (1, 0, . . . , 0) G , and where i?„,_i(0, i5) denotes 
the Euclidean ball of R"^^ of center and radius D. Using spherical coordinates 
z' = \(x - x')^\ei+rn, (r,fi) G (0, +00) x §""2 ^^^^ ^ _ (sin(^), cos(^)e) ((^,6) G 
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(— |, |) X and using the estimate — x')_L|ei + rQ\ > |r + |(x — a;')_L| sin 

we obtain 



2 



< / " ( ] ^ Z I ■ ..M^ ^" 1 (3.72) 

Note that 

"20 (ji pSD 



rdr < / — r(ir < oo, (3.73) 
lo |r + |(x — x')_l| sin(6')| 2 J~3D\r\2 

for ^ G (-f , f ). Combining (1X70]) . flX72D and flX75]) . we obtain 

nil I 3, ^ 3^ 

W,„(x,a;',i;') < ^ ^/ cos"^'(e)rf^ / r-^rfr. (3.74) 

Therefore fl3.55p holds for "m"= m + 1 = n. □ 

Proof of Lemma \3.3[ We first prove the estimates fl3.75p and fl3.76p given below 

/ I <C2-C'2 \n{\zo -z\), when n = 2 (3.75) 

/ ; — ; ; — 7 < ; 7, when n > 3, (3.76) 

Jx\zo- zi\^\zi- z\^-^ - \zo- z\^-^' 

for (2:0, z) G and for m G N such that Zq ^ z and 2 < m < n — 1, where the positive 
constants C'2 do not depend on {zq,z). 

Let (20 ; -2^) £ X^ and let m G N be such that zq ^ z and 1 < m < n — 1. Performing 
the change of variables zi = z + ri^i, (ri,r2i) G (0,/^) x S"^^ (where D denotes the 
diameter of X), we obtain 

<-rf ^.^^n-n .-^^n(3.rT, 



I X \zo - zi\'^\zi - z\^ 1 Jo i§n-i |rfi - ril^il™ Jq Jgn-i |rei - nfin 
where ei = (1, 0, . . . , 0) G M" (we also used a rotation that maps VL to the vector ei) and 

r = l^-^ol, 1^= (3.78) 

Performing the change of variables f2i = (sin(6'i), cos(6'i)6i), (6*1, 0i) G (— |, |) x S""^), 
on the right hand side of (13.771) we obtain 

dn = c{7i) / / — —TTTTZ^dri, (3.79) 




^S" 



1 |rei — ril]i|™ Jq J _2l (r^ + — 2rri sin(6'i)) 2 



where c(n) := (by convention c(2) := 2). 

Consider the case n = 2 and m = 1. Using (13 .79^ and the estimate (r2 + ri2 — 
2rri sin(6'i))^ > 2^^(\ri — rsin(^i)| + r| cos(^i)|), we obtain 



Sir) := r r . ^ rfr, (3.80) 

Jo (^^ + n - 2rri sm(6'i)) 2 
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< r . 



Jo |ri - r sin(6'i)| + r| cos(6'i)| 

< r r < r r^^^^ 

Jo i-2D ril +r|cos(^i)| Jq Jo ri + r|cos(^i)| 







r|cos(6'i)| / Jq \r|cos(6'i)| 



(3.81) 

where Ci := 2t ^(ln(4D) -ln(| 003(^2) |))t^^2 < 00. Estimate (I3T5D follows from flSTTj) . 
fl3T9ll and flaTSTD . 

Consider the case n > 3 and 2 < m < n — 1. Note that 



r cos(ei) < y r2 + - 2rri sin(^i), (3.82) 
for (r,ri,^i) G (0, +cx))2 x (-f, f). Combining fl3T9D and ^M), we obtain 

Jo \rei - nf^il™ - r™"2 (^2 + ^2 _ 3^^^ sin(^i)) ^ 

( 



<nr f -r^''^'"' dr, ^ n r ^^^^^ 



j,m—2 



'0 



2 



+ — 2rri sin(0i)) r™ ^ Jg '"i 



r,m— 1 



rfr/ (3.83) 

7] 

(we perform the change of variables ri = r?7, dri = rdrj). Finally fl3.76p follows from 

In ( 1+'? 

(I3:83|) and fl3T71) and the estimate J^^°^ dr] < +00. 

We are now ready to prove ( I3.14p . Let n > 2. From (13. 8p . it follows that 

i?(2o,2:n+i), (3.84) 

for a.e. [zq, Zn-\-i,Vn+i ) eX X X X §"-\ where 

R{zo, Zn+i) = I —TT. — 7—[dzi ■ ■ ■ dzn, (3.85) 

JX" ffi=l \Zi- Zi^i^' 

for (20, Zn+i) e X X X. 

Assume n = 2. Combining (13.751) and fl3.85p we obtain 

RU^,)< (3.86) 

JX 1^1 ~ ^o| 

Performing the change of variables zi = zq + riQi, on the right hand side of f l3.86p . we 
obtain ^ 

R{zo,Z3)< [ [ {C-C'\n{\r,n, + zo-zs\))dn,dr,. (3.87) 
Jo Js^ 
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Assume n > 3. Using fl3.85l) and fl3.76p . we obtain 

R{zo, Zn+i) <C I -— -— -dzndzn-i, (3.88) 

Jx'^ l^n+1 ~ Zn\ \Zn ~ ^n-l| |^n-l ~ Zq\ 

where C does not depend on {zo,Zn+i)- Performing the change of variables Zi = 
Zi+i + rifli, {ri, Qi) G (0, +00) x §"~^, i = n — l,n, we obtain 

drn-idvndQn-idQn 



(0,D)2x(§"-i)2 \Zn+l — Zq + Vn^n + T"n~l^n~l\ 

drn-idQn-i 



{0,D)xS"-i J{0,D)xS 



571-1 \\zn+i — Zq + r„r2„,|ei + r„„i 



-drndQn, (3.89) 



where ei = (1, 0, . . . , 0) G W^. Performing the change of variables Q = (sin(6'), cos(6')0), 



{9, 9) G (-f , f ) X §"-2^ we obtain 

r f II 01 ^^ ^ 4n)S{\w\), (3.90) 

for U7 G M", w 7^ 0, where c(n) := |S"^^| and 5'(|w|) is defined by (13.801) . Therefore from 
(ICTj) . (IXHTD and (IXHUD . it follows that 

i?(2;0,2;„+l) < / (C - C"ln(|z„+i - Zo + '"nfin|))c?'^nC?fin, (3.91) 

-'(0,D)x§"-i 

where the positive constants C, C do not depend on {zQ,Zn+i). 
Finally from flCTl) and flMB . it follows that 

i?(2;o,2;„+i) < / (C - C"ln(|z„+i - zo + '"n^n|))(ir„(iO„, (3.92) 

-'(0,D)x§"-i 

for n > 2 and for (zq, Zn+i) G X'^, Zq ^ Zn+i, where the positive constants C, C do not 
depend on {zQ,Zn+i). 

Let n > 2 and let {zq, Zn+i) G X^, zq ^ Zn+i- From (13.921) and the estimate 

\rn + (^0 - Zn+l) ■ < k„fi„ + Zq - Zn+l\ it foUoWS that 

R{zo, Zn+l) < c{n)CD - C / ln(|r„ + (^o - z^+i) ■ fi„|)rfr„rffi^ 

^S"-i Jo 

/2D 1-20 
\n{\rn\)drn) = c{n){CD - 2C' / Hrn)drn), (3.93) 
2D Jo 

where c(n) := |§"-i|. Statement fl37[il) follows from and ([333]). 

We prove (I3.15p . We first obtain 

K''^\l-K)-^J(f){x)= [ -fn+i{x,x',v'){I - K)~^J(j){x,v)dv, (3.94) 

for G X and for G L^(r_, c/^). Therefore using (I3.14p we obtain 

_ J^)-i J0(x)|Uoo(x) < C||0|Ui(r_,,5), (3.95) 
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for G L\T^,d^), where C := ||7n+i||L-(XxXx§"-i)||(^ - ^0"^^IU(Li(r_,d€),iH^xS"-i))- 
Therefore there exists a (unique) function $ G L°^(X x r_) such that 

K^'+^I - KY^J(t){x) = j ^x,x',v')(t){x',v')di{x',v'). (3.96) 

Set Vn+i{x,x' ,v') := \i^{x') ■ v'\^{x,x' ,v') for a.e. {x,x',v') E X xT_ and recall the 
definition of d^. Then (I3.15P follows from (]3.96p .D 

4. Derivation of stability estimates 
4.I. Estimates for ballistic part 

Proof of Theorem \EE Let G L°°(X), < 1 and G Li(r_,c/0, 

ll^lUnr-.^O < 1- We have : 









/ 0(x) 




{x)dx 


Jx 







< \\A-A\\c(LHT^40,L-ix)). (4.1) 
Using (14.11) and the decomposition of the albedo operator (see Theorem 12.11) we obtain 

|Ao(0,^)| < \\A-A\\c^L,^r^4i),L^x)) + \^i{(t^M. (4-2) 

where 

Ao(0,V^)= / (l){x){aa{x,v')E{x,x - T-{x,v')v') 

—aa{x, v')E{x, X — T_ (x, v')v')^ ip{x — T_ (x, v')v', v')dv'dx, (4.3) 



Ai(0,V^)= / (f){x) / {Ti-Ti){x,x,v')ip{x,v')dfi{x')dv'dx, (4.4) 
Jx Jv- 

and where Fi = ai + r2 (and Fi = ai + r2). 

Note that performing the change of variables x = x' + tv' on the right-hand side of 
(14. 3p . we obtain 

Ao(0,^)= / \{x' + tv'){r]-fj){t;x',v')dttP{x',v')d^{x',v'), (4.5) 



r_ Jo 



where 

r]{t;x\v') = aa{x' + tv\v')e-^o-i-'+^-'y)ds^ fi{t;x',v') = aa{x' + tv' ,v')e- . 

(4.6) 

We use the following result: For any function G G L^{r^,dC,) and for a.e. 
{xq,v'q) G F_ there exists a sequence of functions ipe^x'^y^ G L^(X-,d^) (which does 
not depend on the function G), \\^e,x'o,v'J\L^r.,d^) = 1, '^e,x'o,v'o ^ ^ supp?/;^,^;^,^^ C 
{(x', v') G F_ I |x' — XqI + |w — WqI < e} such that the following limit holds 

Gix',v')ij,^^>y^ix',v')dax\v') Gix'o,v'o), as £ ^ 0+ (4.7) 
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(we refer the reader to fTH Corollary 4.2] for the proof of this statement). In particular, 
the limits (14.71) holds for some sequence ipe,x[^,v[^ when (xq, v'^) belongs to the Lebesgue 
set of G denoted by C{G) and the complement of C{G) is a negligible subset of r_ (see 

my 

Let C := \^menyj{o}t^{Gm) where Gm is the measurable function on r_ defined by 
G^{x',v') = r^''^\"'{v-v){t;x',v')dt, (4.8) 



for m e NU{0}. Note that the complement of £ is still negligible. Let {xq, iQ E C and let 
(f) G Co(0, t_|_(xq, v'q)) where Co(0, t+(xq, iQ) denotes the set of continuous and compactly 
supported functions on (0, r+(xQ, fg)). Consider the sequence (0m) G (L°°(X))^ defined 
by 

0^(x) = X[o,^)(l^1)0W, (4.9) 

L ' m + l ' 

for X e X and m G N where x = x'n + ivn + x', x' ■ v'q = 0, where X\o ^— 1(^) = 
when t > and X[o,^-^)(^) = 1 otherwise. The support of (pm concentrates around 
the line which passes through x'q with direction v'q. From (14.51) . (14.71) . (14. 8 p (and the 
Stone- Weierstrass theorem), we obtain that 

lim Ao(0m,^£,x' «') = / (f)m{x'Q + tv'Q){r] -fi)(t;x'Q,v'Q)dt. (4.10) 



For the single-scattering part, using (I2.2ip and (13.91) and (13. lip for m = 1 we obtain: 



|Ai(0m,V^e,xj„.^,)|<C ^Ux',v')^,y^,,>{x',v')dax',v') (4.11) 
where C = 11 1 ^J^]^^n^^^f \^ Woo and $ is the function from r_ to M defined by 

^rn{x',v')= / UJn{x,x',v')\(f)m{x)\dx (4.12) 



JX 

(where Wn{x, x', v') is defined in Theorem 12. 4p . From the definition of Wn, it follows that 
$m is a bounded and continuous function on r_. Therefore using Lebesgue dominated 
convergence theorem, we obtain 

lim / ^rn{x\v')'lps,x'v'{x',v')d^{x',v') = ^rnix'o,v'Q)= / W„, (x, Xg , Wq) 1 0„ (x) | rfx. 

£^0+ 7r_ Jx 

(4.13) 

Then using (14.90 . (I4.13P and Lebesgue convergence theorem, we obtain 

limm^+oo fx u!n{x, x'q, v'Q)\(f>m{x)\dx = 0. Therefore taking account of (14.111) and (I4.13p . 

we obtain 

lim lim Ai{^rn,ipe,x'v') =0. (4.14) 

m— >+oo e^0+ 

Combining (g^]) (with = ipe,x',v'), flCTj) . (1^ . (BHD, we obtain 

fT^. (a;',i)') 



(j){t){7] -f]){t;x',v')dt 



< \\A- A\\cimr.,di),LHx)), (4.15) 
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for {x',v') G C and G Co{0,t^{x' ,v')). From the density of Co(0, r+(x', w')) in 
L^(0,r+(x',t/)), it follows that fH?T5D also holds for (j) G L^(0, r+(x', t/)). Applying 
f l4.15p on := sign(r7 — fi){t; x', v') for a.e. t G (0, r+(x', f')) (where sign(s) = 1 when 
s > and sign(s) = —1 otherwise), we obtain (12.231) . □ 



Proof of Theorem 1^.51 Let us define 
C(t) ■.= C{t-x\v') =e-/o-''(^'+««'.^')rf«, r/(t) :=r/(t;x',w') = (Ta{x' + tv' ,v')C{t-x\v'). 

Note that r/(t) = ^C(^)- Then, 



^(C-C)(s)c?s < / \ri-fi\{s)ds<\\A~A\\c[L^i^V-,di);L^{x)). 



thanks to (I2.23p . The point-wise (in t) control on C{t) and the estimate (I2.23P for ri{t) 
show by application of the triangle inequality and the fact that C,{t) is bounded from 
below by the positive constant ^-^^^^'■^'^i^) that 



T+{x' ,v') 



aa{x' + tv\ v') - aa{x' + tv', v') dt <C\\A- i|U(Li(r_,d€);Li(x))- (4-16) 



Here, the constant C depends on M and is of order e^<iiam(x)_ j^^^ surprisingly, 
reconstructions deteriorate when the optical depth Mdiam(X) of the domain increases. 
This shows that 



X 



aix,v') - aaix,v')\dx < C\\A - A\\c(LHr-,dO;LHx))- 



This concludes the proof of the theorem. □ 



Proof of TheoremlKE Let G := {{t, x', v') G (0, +oo) x r_ | t G (0, r+(x', v'))} and let 
h G L°°{g) be defined by 

j-t r-T+{y',v') 

h{t] x' , v') = - a{x' + sv',v')ds + a{y' + sv' ,v')ds, (4.17) 

Jo Jt 

for (t,x',v') G Q. The function h G L°°{Q) is defined similarly. 

We first prove a stability estimate (I4.23P on h, h. Let {yo,vo) G r_, t G 
{0,T+{yo,vo)). Set {yi,vi) := (yo + r+{yo,vo)vo, -vq) and ti = T+{yo,vo) - t. Due 
to the symmetry of cTq, a with respect to the speed variable v (ua, bounded from 
below), we obtain: 

V{tuyi,vi) ' r]{ti;yi,vi) 

where 77 is defined by (14.60 . Note that from (I4.17p . it follows that \h(t;yQ,VQ)\ < 
||(T||ooT4-(yo! "^o) < -D||o"||oo where D denotes the diameter of X. A similar estimate 
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is valid for h and a. Therefore using the fact that \a — a 
a = h{t), we obtain 



\h-h\{t;yo,Vo) < e^'°'^-(li-ll-'ll'^il-) 



^ g- min(a,a) jgO _ ga| 



V{ti;yi,vi) r]{ti;yi,vi) 



(4.19) 



Note that from (I4.6P it follows that 

< aoe-^ll"!l°° < r]{s; y, v) < (4.20) 
for {y,v) G r_ and s G {0,T+{y,v)), where D is the diameter of X. A similar estimate 

hi 

bb 



is valid for (aa, rj). Therefore using the equality f — f = "'"'^'•'^ (for a = ri{t; yo, vq), 



b = rfiti] yi, fi)), we obtain 

ri{t;yo,Vo) v{t;yo,Vo) 



< {v{t;yo,Vo)\v - V\iti;yi,vi) + r]{ti;yi,vi)\r] - r]\{t;yo,Vo) 



< 



,Di\\a\\^ + \\a\\oo) 



WaWoolv - v\i'ti;yi,vi) + lka||oo|?7-?7l(t; 2/0,^0))- (4.21) 



Using (14.211) . integrating in the t variable (ti = T+{yo,vo) — t) and using (12.231) . we 
obtain 



r]{t;yo,Vo) r]{t;yo,Vo) 



V{ti;yi,vi) r]{ti;yi,vi) 



gOdkii^oo+lklloo) 

dt < ( IICT, 



a 1 1 oo ~r II 1 1 oQ 



X \\A-A\\c(LHr.40,LHx))- (4-22) 



Combining (14.191) and (I4.22p we obtain 



\h - h\{t;yo,vo)dt < C\\A - A\\c{L^r-,d^),LHx)) 



(4.23) 



for a.e. {yo,Vo) G r_, where the constant C depends only on D, ||cr||oo, ||5'||oo5 ||c"a||oo 
and 1 1 a"a I loo- We now prove (14.261) . Estimate (I4.26P is, in particular, a consequence of 
the identities 

dh{t;yo,Vo) o / , n dh{t;yo,Vo) . 

= -2a{yo + svo, vq), = -2a{yo + svq, vq), 

at at 

for a.e. {yo,vo) G r„ and t G (0, r+(?/o, ^^o))- 

Consider the operator Tsa{.,v) associated with da = a — a : 



(f) ^ 5a{x,v)(p{x)dx 



X 



for V G S" ^. Using the change of variable x = y + tv {dx = \i^{y) ■ v\dtdy, 
y G dX.{v) := {x' G dX \ u{x') - v <0},t e {0,T+{y,v))) we have 



\{T5a(.,v),V) \ = 


/ W{y)-v\ 




JdX-(v) 



Sa{y + tv, v)(p{y + tv)dtdfi{y) 
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and using integration by part in the inner integral we have: 



\{T5ct{.,v),^)\ < 



Wiy) ■ v\ 



{h — h){t; y, v)v ■ V(p{y + tv)dt 



dfi{y) (4.24) 



< 211^^11^ 



h-h 



(t; y, v)dtdfi{y) 



ldX-(v) Jo 

< C\\A - A\\ciLHr.,dO,LHx))\\<^\\w^'°°{x)^ 
for a.e. v G S"~^ (we used 04.231) ). Identifying 6a and Ts^, we obtain 

Ik — o-||l°o(§"-i,vk-i.i('^)) — ~ •^\\ciL^ir-,d^),L^{X))- 

We prove the following estimate 

Ika — 5'a||L°°(§"-i,Li(X)) < C\\A — A\\ c(L^(T -);L^(X)) , 



(4.25) 



(4.26) 



(4.27) 



where C depends on the uniform bounds cto and M. Combining (14.261) and (14.271) . we 
obtain (I2.28p . which completes the proof of Theorem 12.61 The estimate on aa{x,v') 
comes from an estimate on cr(x' + sv',v')ds by the triangle inequality. The latter is 
given by 



2 / a{x' + sv',v')ds 



T+(x',v') Pt 

a{x' + sv' ,v')ds + i / - 
^Jo Jt 



T+{X',V') 



jcr(a;' + sv', v')ds. 



We just obtained control of 



T+{X',V') 



j ((J — 0") {x + sv' , v')ds, 



in the sense in the t— variable. It thus remains to control the constant term 
£+^'''''''\a-a){x' + sv',v')ds. 

Note that the latter term is nothing but the X-ray transform (Radon transform 
when = 2) of 0" along the line of direction v' passing through x'. In the setting 
of measurements that are supposed to be accurate in C{L^{T^,dC,); L^{X)), the line 
integral is not directly captured as it corresponds to a measurement performed at a 
point X = x' + T^{x' ,v')v' . This is the reason why we assume that cr is known in the 
5— vicinity of dX. 

Knowledge of a and a in the 5o- vicinity of dX allows one to control Jq^^^ ^ cri^x' + 
sv',v')ds by \\A — A\\c{L'^{r_,d^):L^{x))- When the X-ray transform of a is well captured 
by available measurements, as for instance in the presence of boundary measurements 
[2], then Sq can be set to 0. 

More precisely, we find that 

/•T+{x',v') /•T+{x',v') 

/ {a — d'){x' + sv' ,v')ds = / (f){s){a — d'){x' + sv' ,v')ds, 

Jo Jo 
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where G C^{0,t+{x',v')) is equal to 1 ior Sq < s < r+(x',f') — Sq. Integrating by 
parts, this shows that 

T-\-{x' ,v') 

(cr — cr)(x' + sv', v')ds 
< \\{a - + tw',w')IU-i.i(o,r+(x'y))ll'^ll^''^(0'^+(^''^'')) ^'^ll^ - M\c(L\v.M);LHx)). 
thanks to estimate fl4.23p . By Lipschitz regularity of the exponential, we thus have that 

rr+{x',v') 

Jo 



dt<C\\A-A\\ciL^ 



c(L^T-,dO;LHx))- 



(4.28) 



The stability result fl4.27p on o"a(x, v') follows from (12.231) and the triangle inequality as 
in the proof of Theorem 12. 5[ □ 



Proof of Corollary \2. 7[ Let p > 1. We first derive (12.291) from the following estimate 



|cr — cr||ioo(gn-i^^y-i,p(jf)) < C\\A — AYnix 



(4.29) 



Assuming that a - a is bounded by Cq in (§"-\ W^p{X)), then using (ICTj) and 
using the complex interpolation result [5l Theorem 6.4.5 pp. 153], we obtain 



\a — a\ 



<c'\\A-A\\r';r. 



C(L^(T_,d£.).LHX)y 



(4.30) 



for -1 < s < r, where s = (1 - ^) x (-1) + rO and C = C^-^C^ This proves flZ^ . 

We prove (ICTD . From KIT} , it follows that \h-h\{t; x' ,v') < {\h\ + \h\){t;y o,Vo) < 
-D(||cr||oo + llcrlloo), where D is the diameter of X. Hence we obtain 



h-h 



{t;x',v') < DP-\\\a\\^ + ||^|U)p-> - h\{t;x',v'), 



(4.31) 



for {x',v') G r_ and t G (0, r+(x', f')). Using (I4.24p and Holder inequality, and using 
(I4.3ip . we obtain 



(cr — cr)(x, v)ip{x)dx 



X 



1 

< - 

- 2 



T+ (y,v) 



\h-h\P{t;y,v)dtdM HVvJlLp'm 

dX.{v)JO J 



< C\\-^ - ■^\\c{L^r-,d^),LHx))\\'^\\w^'p'(x)y (4.32) 

for a.e. v G S"-^ and for G (X) := {tP G L^' (X) \ supp^/; C X, V^A G Lp'{X, C")} 
where p'~^ + = 1. Estimate (14321) proves (l4l29|) .D 



4-2. Estimates for single scattering 



Proof of Theorem \2.4\ Let (x, x') E X x dX. Set v' = x — x' and let v' G S" ^ be 
such that v' ■ w'"*" = 0. Let = |x — then x = x' + t'^v'. First assume n = 2. Using 
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ai - 


- «! (x + ef'"*", x', v') 




i^(a;') ■ v' 


(l + ln 


(\(t'„-T+(x',v'))v'+ev'^\+iT+[x',v'}-t'„)\ 


) 


[ \t'ov'+ev'^\-t[, ) 



< 



{T,-f,){z,z',w') 



|z/(z') ■ W'\W2{Z, Z\ W') 



1 




oo 


\u{z>)-w'\ 




{\{t[,-T+(x',v'))v'+ev'^\ + {r+(x',v')-t'^)\ 


y \t'^v'+ev'^\-t'^ ) 



(4.33) 



for £ > 0. Therefore using fICT]) as e 0+ and flT^ . we obtain ^TI^. 

Assume n > 3. Using the equahty cti = Fi — and using (12.211) . (13.1 ip (for 
"m = 1"), we obtain 



|«i - tti I (a; + ev'^, x\ v') ^ 



+ e 



{T,-f{){z,z',w') 



L°°(xxr_) 



\z — z' — ((z — z') ■ 


w')w' 


''-H^2-T2){Z,Z',W') 




u{z') ■ 


w' 





L°°{xxr_) 

for e>0. Therefore using as e ^ 0+ and fl2:22|) . we obtain fl2:25|) . 



(4.34) 
□ 



Proof of Lemma \2. 9\ . First consider the case n = 2. The statement in Lemma [2.91 is a 
straightforward consequence of (12.351) and we have the following inversion formula: 



(4.35) 



Let now n = 3. Using (12.361) and the identity cos(6') = 1 — 2sin^(|), 6* G M, we obtain 



^9 



sm 



(o) dO^ 



for g G [0, 1). Performing the change of variables 9 = 2 arcsin(;^7=y) {d6 
obtain 



2dt 



7j, we 



27rhig) 



1+g 



(l + t2)l( 1 + 



Performing the change of variables v = jz^t, we obtain 



i + g 

l~9 



e dt. 



2Tih{g) 



+ 00 



I 



2 / \ 2 
1+g ' 



\ (l+t;2) 



dv 



Y r+ca 



dv, (4.36) 



for g G [0, 1). Note that from the above, it follows that 27ih{g) = + o{j^) a.s g ^ V 



(where we used the integral value f^"^ — 

Differentiating (I4.36P with respect to g, we obtain 
1 /■+°° dv 



+ 00 



!)• 



2'Kh{g) 



(l-^7)'7o il + v^)'2y/ujj{g)+ujl{g)v^ {I + 9} 



rM9), 



(4.37) 
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for g e [0, 1), where h = ^, uji{g) := U2ig) := j^, and 



h2{g) ■■- 

Integrating by parts, we obtain 

P+OO 

h2{9) 



+00 



v'^dv 



{l + v^)-^^ul{g)+^l{g)v^' 
ujl{g)dv 



(4.38) 



(4.39) 



{1 + v^)Hu;l{g) + ul{g)v^)-- 

for g G [0,1), where we use that a primitive of the function r(f) := v{l + w^)~t is 
given by the function R{v) := —(1 + v'^)~i and where we used that the derivative of 

1 3 

the function s{v) := {ujl{g) + u;|((yf)t>^) ^ is given by s{v) = uj^g) {ujl{g) + ci;|((7)t>^) ^ . 
Performing the change of variables "f"= j^v on the right hand side of (14.391) . we obtain 



h{g)=ujl{g)£^{l + v'rl fl + 



dv. 



(4.40) 



for g e [0, 1). Combining (14.371) and (14.401) . we obtain 



2nh{g) 



(l + t;')"2 1 + 



^-9 



v'^\ ^ dv 



1-9 
1+g 



+ CXJ 



{i+v'r-2(i + 



1+9 



' dv 



(4.41) 



for g e [0, 1). Using the estimate 1 + (izf ) > 1 + (j^^ for v G (0, +oo) and 
g G (0, 1), we obtain that the second integral on the right hand side of (14.411) is less than 
the first integral on the right-hand side of (I4.4ip for g G (0, 1). Therefore using also 

that the second integral is multiplied by (^i^j (< 1); obtain h{g) > for G (0, 1), 
which proves that h is strictly increasing on (0, 1). □ 

Proof of Theorem {MUU We prove fl237l) . From I^M), dZSS]), fl230D and fl23TD . it 
follows that 



\E{x,x')ag{x) — E{x,x')ag{x)\ < 



{T^-f,)ix,x',v') 



\iy{x') ■ v'\uJn{x, x', v') 



(4.42) 



for {x,x') G X X r_. Using (14.421) and the estimate mm{E{x,x'), E{x,x')) > 



-Dmax(||o-||oo,||5-||oo) 



'see (12.181) ) where D is the diameter of X, we obtain 



(ri-fi)(a;,a;',^') 



|z/(a;') ■ v'\wn{x, x', v') 



+ \E-E\{x,x')), (4.43) 



for (x, x') G X X r_ (we used the identity ab — db = {a — d)b+ {b — b)d for a := Cg^x) and 
b := E{x,x')). Let {x',v') G r_. Integrating (I4.43P over the line which passes through 
x' with direction v' and using the stability estimate (I4.28P and (I2.18p . we obtain 

(ri-fi)(x,x',t;') 



T+{x' ,v') 

\ag - ag\{x' + tv')dt < C 



\l'{x') ■ v'\Wn{x, X', V') 



+ - A|U(ii(r_,rf5),Li(x)),(4.44) 
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where C = 6^™^^^""^!!°°'"°"!!°°) max(i5, C) and C is the constant on the right-hand side 
of fl4.28p . Integrating fl4.44p over dX_{v') := {x' E dX \ v' ■ u{x') < 0} with measure 
lu^x') ■ v'\d^{x'), we obtain (12.371) . We now prove (12.381) . Assume ||/i((7)||t^i,oo(x) < 
ll^(^)l|iyi'°°(x)- Using (12.341) and min(crs,5"s) > Cs^, we obtain 

\\h{g) - h{g)\\Liix) < —(\\i(^s - &s)h{g)\\mx) + \Wg - agWmx)) (4.45) 

(we used the identity ah — ah = [a — d)h + [h — h)a for a = as and h = h{g)). Using the 
identity cTg = o" — (Xa, we have 

\\{as - &s)h{g)\\Li{x) < \Wa - ^a\\mx)\\h{g)\\L^{x) + \\{(^- &)h{g)\\LHx) 

< IWa - ^a\\L^X)\\h{g)\\Lo-(X) + IW - cr || || V/l(5') || L^CX) (4.46) 

(we used the fact that a = a at the vicinity of the boundary dX). Finally combining 
dHSD, fl06l) . (EJZD and (E^H]), we obtain fl238l) .D 

Appendix A. 

For m > 2, let (3m denotes the distributional kernel of the operator where K is 
defined by (13.21) . We first give the explicit expression of P2, Ps- Then by induction we 
give the explicit expression of the kernel f3rn- Finally we prove Lemma [3.1[ 
From (13. 2p it follows that 

PT—(x,v) P PT-{x — tV,Vl) 

K'^'ip{x,v) = / / k{x — tv,vi,v) / E{x, X — tv, X — tv — tiVi) 



Jo J§"-i Jo 

X / k{x — tv — tiVi,v' ,vi)iIj{x — tv — tiVi,v')dv'dtidvidt, (A.l) 



,1-1 



for a.e. {x, v) E X xS^ ^ and for ip G L^{X x §" ^). Performing the change of variables 



■ tv — tiVi {dx' = ti ^dtidvi) on the right hand side of ( lA.ip . we obtain 

K'^'ijj{x,v) = / / k{x — tv, X — tv — x' ,v)E{x, X — tv, x') 



Jx 

X / k{x', v', X - tv - x')xIj{x\ v')dv'dx'dt, (A. 2) 

for a.e. {x,v) eX x §"-1 and for e L^{X x §""!). Therefore 

r'<'~{x,'") k{x — tv, Vi, v)k{x', v', Vi). ^^z::^^, 
/32{x,v,x',v') = E{x,x-tv,x) ■ " dt, (A.3) 

for a.e. {x,v,x',v') eX x §"-i x X x §""1. From ([Ql) and it follows that 

K^tjj{x,v) = / E{x,x — tv)k{x — tv,Vi,v) / (32{x — tv,vi,x' ,v')il){x' ,v')dx' dv' dvidt 

i(0,T_(x,i)))x§"-i ixx§"-i 

r-(^'^U r-^^-'"^^^^^ E{x,x-tv,x-tv-t,v,,x') 

ip[x ,v ) / / k[x — tv,vi,v) 



X = X 



f,l I n — 1 



X /i;(x - tt; - txvx, V2, Vi)k{x', v', V2)^^^^^^^'::^^_^,dtidvidtdx' dv' , (A.4) 
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for a.e. {x,v) G X x and for ip G L^{X x S"~^). Therefore performing the change 
of variables z = x — tv — tiVi {dz = t^~^dtidvi) on the right hand side of (1A.3I) we obtain 

■.T-{x,v) n E{x,x -tv,z,x')k{x',v',z - x') 



/33{x,v,x',v' 



Jx 



\x — tv — z\^~^\z — x'\"-~^ 



X k{x — tv,x — tv — z, v)k{z, z — x',x — tv — z)dtdz, (A. 5) 

for a.e. {x, v, x', f ') G X x S"^-'^ x X x S""^. Then by induction we have 

^ . ^ r-^"'"^ f E{x,X-tV,Z2,...,Zm) 

Jo JX^~^ \x-tv-Z2r' \Zi-Zi+i\^ 1 

X k{x - tv, Vi, v)UZ2HZi, Vi, .,=..^1, i=2...m-ldtdZ2 • • • dZm-1, 

for a.e. {x,v,x',v') G X x S^^i x X x §"-^ 

We prove ([X7D for m = 2. From fOOj) . (1^ and ([S3]) it follows that 

/• /• fr-{x,v) k(x — tv,vuv)k(y,v',v-i), --r~~ 

J§"-i Jxx§"-i 7o |a; — tf — y| 

X E{y, y - T_{y, v')v')ip{y - T_{y, v')v', v')dydv'dv, (A.7) 

for a.e. x G X and for ip G L^(r_, (i^). Then performing the change of variable z = x—tv 
{t = \x — z\, V = X — z and dz = f^'^dtdv) on the right hand side of (1A.7I) . we obtain 



K'^Jip{x)= / aa{x,x — z)E{x, z,y) 



k{z,vi,v)k{y,v',vi)^^^^-^ 

, , I ry* 'v I TL 1 I 'v _____ It I I 1 

'xxxx§"-i F z\ \z y\ 

X ^(y, y - r_ (y, - ^- {y^ v')v\ v')dzdydv', (A.8) 

for a.e. x G X and for ■?/; G L^(r_, c?^). Performing the change of variables y = x' + t'v' 
(x' G 9X, t' > 0, (i2 = |z^(x') ■ on the right hand side of (lA.Sp . we obtain 

K'^Jip{x) = Jp a2{x,x' ,v')ip{x' ,v')dfi{x')dv' for a.e. x G X and for ip G L-^(r_,(i^), 
which proves (13. 7p for m = 2. 



Then we prove (13. 8p before proving (13. 7p for m > 3. Let m > 3. From (IA.6P and 
the definition of the operator K'"^ (13. 6p . it follows that 

K'^lpiZo) = (ra{Zo,Vo) f]m{Zo,Vo,Zm,Vm)i^{Zm,Vm)dz^dv^dVo 

(^aiZo^VQ) II I E{zo,Zo-tVo,Z2,...,Zm) 

X 



X dt(i2;2 . . . dz^_i4j{zm, Vm)dzmdvmdvo, (A. 9) 

for a.e. zq & X and for G L^(X x §"~^). Therefore performing the change of variables 
zi = zq — tvf) {dzi = f^'^dtdvo, t = \zo — zi\ and vq = zq — zi) on the right hand side of 
dXH) we obtain K""iIj{zo) = Jj^^^r.-i 7m{zo, Zm,Vm)'4'izm,Vm) dzmdvm, for a.e. 2:0 G X 
and for tp e L^{X x S"-"^), which proves fl3l8|) . 
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We prove ([X7D- Let m > 3. From ([SSD and it follows that 



_ f f E{Zo, . . .,Zm,Zm- T_{Zm,Vm)Vm) , , _ ^ ( „ „ WAin^l 

X [o-a(2;0, ^'o)n™ iA;(Zi, fj, Vi-.i)\^_^r^r;~^^^ i=0...m-l • • • dZm-ldZmdVm, 

for a.e. ^ and for G L^(r„,(i,^). Then performing the change of variables 

Zm = "^m" +t'Vm ( "2™" G t' > 0, dz^ = ■ Vm\dflC'z^'' )dt') , We 

obtain K"^Jip{zo) = J^_am{zo, Zm,Vm)ip{zm,Vm)dfi{zm)dVm ioT a.e. zo e X and for 
e Li(r_,ciO> which proves flXTjl .D 
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